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In a recent paper [Phys. Rev. E 61, 2831 (2000)], Borstnik, Stark and Zumer have studied the 
stability of a colloidal dispersion of micron-sized spherical particles in the isotropic phase of a 
nematogenic material. Close to the nematic transition, the attraction due to a surface-induced 
paranematic order can yield flocculation. Their calculation of the nematic-mediated interaction 
was based on an ansatz for the order-parameter profile. We compare it with an exact numerical 
calculation, showing that their results are qualitatively correct. Besides, we point out that in the 
considered regime, the exact interaction is extremely well approximated by a simple analytical 
formula which is asymptotically exact. 

PACS numbers: 82.70.Dd, 61.30.Cz 




In recent years, a large interest has been devoted to 
understanding the interactions and phase behavior of col- 
loidal particles dispersed in a nematic phase jl], 0, 0, 0, 3] 
or in the isotropic phase of a nematogenic compound 1 1, 
0, ^, 0j. In the nematic phase, colloids experience a spe- 
cific elastic interaction because they induce competing 
distortions of the nematic director field. New physics 
arises due to the long-range character of this interaction 
and the induction of topological defects [[j] . 

In the isotropic phase, the surface of colloidal parti- 
cles can induce a local paranematic order 110 , 
ing rise to a short-range elastic interaction 
effects compete: an attraction due to the favorable over- 
lapping of the paranematic halos and a repulsion due to 
the distortion of the director field. For small particles, 
of size comparable to the nematic-isotropic coherence 
length £, it has been predicted that repulsion may dom- 
inate and stabilize the colloidal dispersion ^|. (Note 
that latex particles as small as 50 nm have been success- 
fully dispersed in lyotropic nematics Q.) On the other 
hand, Borstnik, Stark and Zumer have predicted that for 
micron-sized particles attraction dominates [Q , which al- 
lows to trigger flocculation close to the nematic transi- 
tion g. 

The results of Borstnik, Stark and Zumer S are based 
on a composite ansatz for the nematic director field n 
and for the scalar order-parameter Q, within an uniaxial 
hypothesis. It turns out that our exact calculations |5), 
based on a multipolar expansion for the full tensorial 
order-parameter Qij, rest on the same theoretical model, 
and can be performed also for micron-sized particles. 
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FIG. 1: Paranematic interaction energy between two spher- 
ical particles of radius R — 0.25 /jm as a function of their 
distance to contact d. The parameters are a — 1.8 x 
10 5 Jm- 3 K-\ T* = 313.5K, AT = 1.3K, L x = 9 x 
10" 12 Jnr\ Qs = 0.3, Gq = G n = 10~ 3 Jm~ 2 . The cor- 
responding nematic coherence length is £ = 10.7 nm. The 
dashed-dotted line is extracted from Fig. 7 of Ref. The 
full- line is the exact result, numerically calculated according 
to Ref. j^]. The dotted-line, practically coincident with the 
full-line, corresponds to the asymptotic formula (H). 



In this Brief Report, we check the correctness of the 
paranematic interaction energy used in Ref. |8) , by com- 
paring it with the exact one, numerically calculated ac- 
cording to the method of Ref. || . For the typical values 
considered in Ref. ||, we find that the exact interaction 
is attractive instead of repulsive at distances > 5£, and 
that it is about two times weaker at distances of the or- 
der of £ (Fig. [j]). However, the order of magnitude and 
the sign of the interaction is correct in the regime where 
the paranematic attraction competes with the electro- 
static repulsion. In particular, at particles' separations of 
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the order of £, the paranematic attraction remains much 
larger than the van der Waals attraction: this implies 
that the conclusions of Ref. || regarding critical floccu- 
lation phenomena remain qualitatively correct. Finally, 
we show that in the whole regime of interest, the exact 
interaction calculated numerically is extremely well ap- 
proximated by a simple asymptotic formula [see Eq. (pi)]. 



The Landau-de Gennes 
density considered in Refs. 



paranematic free-energy 
is 



^(VQf + ^xQ^Vnl 2 (1) 



for the bulk, and 

f s = G Q (Q - 
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for the surface. Here Q is the scalar order-parameter and 
n the nematic director. The parameter Q s is the order- 
parameter favored by the particles' surface, and 9 is the 
angle between the direction of n at the surface and the 
normal to the surface v , which is assumed to correspond 
to the easy-axis. The coefficients a and L\ are material 
parameters, AT = T — T* is the difference between the 
actual temperature T and the limit of stability T* of the 
isotropic phase, Gq and G n are introduced to describe 
the strength of the surface anchoring. 

Incidentally, we note that the surface free-energy den- 
sity (0) is inconsistent within a Landau-de Gennes frame- 
work, unless Gq = G n . Indeed, at quadratic order, the 
most general expansion in terms of the tensorial order 
parameter Qij and of the normal v to the surface can be 
written as 



with Q s = — <7i/3<72i- Therefore, Eq. (||) is compatible 
with this expression only if Gq = G n . 

With the required condition Gq — G n = G, the free- 
energy considered in Ref. [|| is identical to the one we 
used in Ref. ||, with the correspondence 
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where we have dagged the quantities appearing in 
Ref. ||. With the above relationships, we have numer- 
ically recalculated the exact interaction energy between 
two spherical particles of radius R as a function of their 
distance to contact d — using the same parameters as in 
Ref. §. The comparison with the results of the ansatz 
of Borstnik, Stark and Zumer M is shown in Fig. |]. 
We find a qualitative agreement, as previously discussed. 
Note that the exact nematic director profile displays a 
Saturn-ring defect P], which is absent in the ansatz of 
Ref. JjJ . For the case of micron-sized particles considered 
in Refs. ^ M, this defect lies however in a region where 
the paranematic is almost completely melted. 

Finally, we have compared our numerical calculation 
with the asymptotic formula obtained by us in Ref. ||. 
With the correspondence (0), the latter reads 



fs — 9i QijViVj + g 2 i QijQij + .922 Q./Qji. l ',i'>- 

+ 32-iQijQkiViVjV k v u (3) 

where summation over repeated indices is implied. For a 
uniaxial tensorial order-parameter 
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where d = d/£, R = R/£, and A is a constant given by 



this yields 



fs = 9i Q + g [3 521 + 2 ( 5 22 + .923 )] Q 2 

- j [2g x Q + (.922 + 4 523 ) Q 2 ] sin 2 9 
+ ^ 5 23 Q 2 sin 4 9. 



(5) 



Matching Eqs. (|) and (H) requires setting 1723 = and 
consequently c/22 = 0. Equation ^ then becomes 
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where I = Li/G^ is the reduced extrapolation length 
of the anchoring. The nematic coherence length is ^ — 
y/3Li/a AT. As shown in Fig. []], the agreement between 
the numerical and the analytical calculations is excel- 
lent in the range of separations relevant to the colloidal 
flocculation discussed in Ref. Q. Owing to its simplic- 
ity and validity, formula (||) thus offers a straightforward 
means to systematically investigate the stability of such 
paranematic-wetted colloids. 



3 



[1] P. Poulin, H. Stark, T. C. Lubensky, and D. A. Weitz, 

Science 275, 1770 (1997). 
[2] P. Poulin, V. Cabuil, and D. A. Weitz, Phys. Rev. Lett. 

79, 4862 (1997). 
[3] H. Stark, J. Stelzer, and R. Bernhard, Eur. Phys. J. B 

10, 515 (1999). 
[4] R. Yamamoto, Phys. Rev. Lett. 87, 075502 (2001). 
[5] D. Andrienko, G. Germano, and M. P. Allen, Phys. Rev. 

E 63, 041701 (2001). 
[6] P. Galatola and J.-B. Fournier, Mol. Crys. Liq. Cryst. 

330, 535 (1999). 
[7] A. Borstnik, H. Stark, and S. Zumer, Phys. Rev. E 60, 

4210 (1999). 



[8] A. Borstnik, H. Stark, and S. Zumer, Phys. Rev. E 61, 
2831 (2000). 

[9] P. Galatola and J.-B. Fournier, Phys. Rev. Lett. 86, 3915 
(2001). 

[10] P. Sheng, Phys. Rev. Lett. 37, 1059 (1976). 

[11] K. Miyano, Phys. Rev. Lett. 43, 51 (1979). 

[12] O. Mondain-Monval, J. C. Dedieu, T. Gulik-Krzywicki, 
and P. Poulin, Eur. Phys. J. B 12, 167 (1999). 

[13] P.-G. de Gennes and J. Prost, The Physics of Liquid 
Crystals (Clarendon, Oxford, 1993). 

[14] Actually, in our free-energy we keep the full biaxial order- 
parameter Qij instead of assuming the uniaxial order (mj . 



